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ABSTRACT 

Aims. We investigate the electron temperature gradient in the galactic disk as measured by young HII regions on the 
basis of radio recombination lines and the corresponding gradient in planetary nebulae (PN) based on [OIII] electron 
temperatures. The main goal is to investigate the time evolution of the electron temperature gradient and of the radial 
abundance gradient, which is essentially a mirror image of the temperature gradient. 

Methods. The recently derived electron temperature gradient from radio recombination lines in HII regions is compared 
with a new determination of the corresponding gradient from planetary nebulae for which the progenitor star ages have 
been determined. 

Results. The newly derived electron temperature gradient for PN with progenitor stars with ages in the 4-5 Gyr range is 
much steeper than the corresponding gradient for HII regions. These electron temperature gradients are converted into 
O/H gradients in order to make comparisons with previous estimates of the flattening rate of the abundance gradient. 
Conclusions. It is concluded that the O/H gradient has flattened out in the past 5 Gyr at a rate of about 0.0094 dex 
kpc -1 Gyr -1 , in good agreement with our previous estimates. 
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1. Introduction 

Radial abundance gradients in the galactic disk and their 
time variations are among the main constraints of chem- 
ical evolution models for the Milky Way. These gradi- 
ents can be determined from a variety of objects, such as 
HII regions, cepheid variables, open clusters and planetary 
nebulae (PN). In a recent series of papers, Maciel et al. 
(|2003l |2TjO"51 [2006]) estimated the time variation of the ra- 
dial abundance gradients taking into account a large sam- 
ple of PN for which abundances of O/H, S/H, Ne/H and 
Ar/H have been derived. Based on individual estimates of 
the progenitor star ages, it was concluded that the radial 
gradients are flattening out at an average rate of about 
0.005 - 0.010 dexkpc" r Gyr _1 for the last 8 Gyr, approxi- 
mately. A comparison of the PN gradients with results from 
HII regions, OB stars and associations, cepheids and, espe- 
cially, open cluster stars, strongly supports these conclu- 
sions. 

On the other hand, it has long been known that a posi- 
tive electron temperature gradient of about 250— 450 K/kpc 
is observed in the galactic disk, mainly on the basis of ra- 
dio recombination line work on HII regions (see for exam- 
ple Churchw ell fc Walmsley [ WSl C hurchwell et al. 
Shaver et al. fTMSl W ink et al. ITTJ8"5l Afflerbach et al. [1591)1 
and Deharveng et al. I2000)) . Such a gradient is interpreted 
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as a reflection of the radial abundance gradient of elements 
such as O/H, S/H, etc. in the galactic disk, since these 
elements are effective coolants of the ionized gas (see for 
example Shaver et al. [1983 ) . 

Recently, Quireza et al. (2006) presented a detailed 
study of a large sample containing over a hundred HII 
regions spanning about 17 kpc in galactocentric distances 
for which accurate electron temperatures were determined 
from radio recombination lines, specifically H91a and 
He91a. The observations were made with the 140 Foot 
telescope of the National Radio Astronomy Observatory 
(NRAO), and are of unprecedented sensitivity compared 
with previous studies. According to this work, the best es- 
timate of the gradient, obtained from a sample of 76 sources 
with high quality data, is dT e /dR ~ 287 ± 46 K/kpc, with 
no significant variations along the galactocentric distances. 
A slightly larger gradient (up to 17%) was obtained by ex- 
cluding some HII regions which are closer to the galactic 
centre, and may not belong to the disk population. 

Regarding planetary nebulae, our earlier work (Maciel 
& Faiindez-Abans 1985} based on a sample of PN classi- 
fied according to the Peimbert types (cf. Peimbert 1978) 
suggested a positive electron temperature gradient in the 
range 550 — 800 K/kpc, somewhat steeper than the HII re- 
gion gradients observed at the time. 

In this work, we take into account the recent PN sam- 
ples analyzed by Maciel et al. (|2003l [20051 [2006)) and derive 
the PN electron temperature gradient for a sample of ob- 
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Fig. 1. Age distribution of the PN progenitor stars in the 
Basic Sample of Maciel et al. ([2006]) . 



jects having similar ages. A comparison of the obtained T e 
gradient with the recently derived value by Quireza et al. 
(2006) for HII regions gives then an independent estimate 
of the time variation of the radial abundance gradients in 
the galactic disk. 



2. The electron temperature gradient in the 
galactic disk 

The determination of abundance gradients is a difficult 
task, basically for three main reasons. First, the magnitudes 
of the gradients are small, amounting at most to a few hun- 
dredths in units of dex/kpc, so that a relatively large galac- 
tocentric baseline is needed in order to obtain meaningful 
results. Second, the uncertainties both in the abundances 
and in the distances contribute to the observed scattering, 
so that large samples are usually needed. Third, chemical 
evolution models generally predict some time variation of 
the gradients, so that it is extremely important to take into 
account in a given sample only objects with similar ages. 
For these reasons, some of the analyses of gradients in the 
literature produce relatively flat gradients (see for example 
Perinotto & Morbidelli I2006)) . On the other hand, accu- 
rate and homogeneous abundances eliminate some of these 
problems, so that relatively steeper gradients are obtained, 
as in Pottasch & Bernard-Salas (20P6"]). 

In our recent work, we made an attempt to overcome 
some of these problems, and estimated the individual ages 
of the PN progenitor stars using an age-metallicity rela- 
tionship which also depends on the galactocentric distance. 
As a result, we have obtained the age distribution shown 
in Fig. [TJ adopting our Basic Sample, which is the largest 
and most complete sample we have considered, containing 
234 nebulae (see Maciel et al. l2"UU3l IMJ51 12H061 for details). 

It can be seen that the ages are strongly peaked around 
4-5 Gyr, where we have 99 objects. In the present work, 
we will consider the objects in this age bracket, in order 
to make comparisons with the younger HII regions. We 
have then collected the electron temperatures of the plane- 
tary nebulae, selecting only the [OIII] temperatures in or- 
der to keep our sample as homogenous as possible. These 
temperatures arc determined from the ratio of the [OIII] 
4363/5007A lines, which are usually among the brightest 
collisionally excited emission lines in the spectra of plane- 
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Fig. 2. Galactocentric variation of the [OIII] electron tem- 
peratures for PN with progenitor ages of 4-5 Gyr. The 
empty circles show some nebulae having extremely hot cen- 
tral stars, not included in the linear regression analysis. 



tary nebulae. We have preferred our own data where avail- 
able (Costa et al. I2004[ I19961 see a list of references in 
Maciel et al . RZOOBI |2"U0"5I [2TO)1 . with ad dition al data by 
Henry et al. @004), Kingsburgh & Barlow (fTSMj) . and Cahn 
et al. (|1992p . The resulting T e variation with galactocentric 
distance R is shown in Fig. [21 where we adopted the same 
distances and solar galactocentric radius as in our previous 
work. The total number of objects in Fig. [2] is somewhat 
lower than shown in the 4-5 Gyr bracket of Fig. [TJ as for a 
few nebulae we could not obtain accurate electron temper- 
atures. 

It can be seen from Fig. [2] that there is a clear tendency 
in the sense that higher electron temperatures are asso- 
ciated with larger galactocentric distances. The best de- 
rived T e gradient for this sample of PN is dT e /dR ~ 670 ± 
65K/kpc, with a correlation coefficient of r ~ 0.76, which 
is similar to the gradient for the 'selected sample' of our 
earlier paper (Maciel & Faundez-Abans I1985p . Adopting 
instead a homogeneous set of [OIII] electron temperatures 
from Henry et al. (2004), which is the largest homogeneous 
sample available for these nebulae, we obtain essentially 
the same result, namely dT e /dR = 680 ± 140K/kpc and 
r ~ 0.65, so that the correlation is real. Our best derived 
slope is illustrated by the dashed line in Fig. [2] 

The average uncertainty in the determination of the 
electron temperatures is generally considered to be within 
10% for the brightest nebulae, which corresponds roughly 
to 1000 K for most objects (see for example Kingsburgh & 
Barlow [SMI and Krabbe & Copetti l2"00l)]) . For the galac- 
tocentric distances an average uncertainty is more difficult 
to establish, as it depends on the adopted distances. Since 
most objects in Fig. [5] are located within about 3 kpc from 
the solar galactocentric radius, an average error of 50% in 
the distances would correspond to a shift in the galacto- 
centric distances of about 0.03 - 1.0 kpc depending on the 
distance and the direction of the line of sight to the neb- 
ula. As a comparison, for the HII regions in the sample by 
Quireza et al. (|2006|) . average formal uncertainties in the 
electron temperatures are within 2%, but systematic errors 
may increase this uncertainty up to 10% for the best data. 
For spectrophotometric distances, average errors of 15% are 
quoted, while for kinematic distances non-circular stream- 
ing motions may increase this figure to about 25%. From 
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the T e x R plot by Quireza et al. (|2006p . an average disper- 
sion of about 2200 K can be obtained, which is about half 
the dispersion in Fig. [2] 

It should be noted that the dispersion observed in Fig. [2] 
is probably real, since the electron temperatures may be 
affected by several factors, such as differences in the ef- 
fective temperature of the central stars, presence of dust, 
optical depth effects, electron density and temperature fluc- 
tuations, etc., apart from the main cause of the T e x R vari- 
ation, namely, the radial abundance gradient. These effects 
are also partially responsible for the observed dispersion 
in HII regions, but for PN the variations in the effective 
temperatures of the central stars and the uncertainties in 
the distances are larger, so that the observed dispersion in 
Fig. [2] is larger than in the case of HII regions. As a con- 
sequence, some nebulae appear not to follow the observed 
correlation very closely. In particular, PN having extremely 
hot central stars, with temperatures in excess of 10 5 K, gen- 
erally have higher electron temperatures than expected by 
their galactocentric distances. Furthermore, these objects 
come from more massive progenitors than most nebulae, 
so that their ages may be lower than 4-5 Gyr, as assumed. 
Some examples include Ml-57, Me2-1, PB6, NGC 6620 and 
a few others, which are plotted in the figure as empty cir- 
cles. Other objects with hot central stars, such as NGC 
2899, NGC 6302, NGC 6537 and NGC 7008 are not plot- 
ted in Fig. [2l as their electron temperatures are too high 
to fit the scale. All these nebulae have not been included in 
the linear regression, so that our derived electron tempera- 
ture gradient applies to stars with temperatures lower than 
10 5 K. In this analysis we have used Zanstra temperatures 
and energy-balance temperatures (see for example Preite- 
Martinez et al. ITMIl Mendez et al. fl"992l Zhang H9931 and 
Stasihska et al. I1997|) . 

Another interesting object is Ml-9, which is the neb- 
ula with the largest galactocentric distance in the sample 
(R ~ 12.4 kpc). In view of its position on the T e x R plane, 
it may single-handedly affect the derived slope. For this 
object, our own results suggest an electron temperature of 
T e ~ 11000 K (Costa et al.[T996), but a more detailed study 
by Tamura & Shibata (fTM?]) and Shibata & Tamura (fT355]) 
gives a larger value, T e ~ 14800 K, which is adopted here. 
This object may then alter the derived slope by about 50 
K/kpc, but the main conclusions of this paper are unaf- 
fected. 

The association of higher electron temperatures with 
lower metallicities can be seen from Fig. [31 where we show 
the inverse correlation between the [OIII] electron temper- 
atures and the O/H abundances for the objects with ages 
in the 4-5 Gyr bracket. PN with central stars hotter than 
10 5 K are also shown as empty circles. Again a relatively 
large dispersion is observed, but the inverse correlation is 
clear, confirming that oxygen is among the main coolants 
in the photoionized gas within the planetary nebulae. 

3. Comparison of the PN and HII gradients 

From a straigthforward comparison of the electron tem- 
perature gradient for PN and HII regions we can already 
conclude that there is some flattening of the gradients dur- 
ing the time interval of about 5 Gyr, which is essentially 
the difference between the ages of the PN progenitor stars 
considered in this work and the much younger HII regions 
in the sample by Quireza et al. (2005]). In other words, the 
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Fig. 3. The inverse correlation between the [OIII] electron 
temperatures and the O/H abundances for PN in the 4- 
5 Gyr age bracket. The empty circles show some nebulae 
having extremely hot central stars. 



conclusions of our recent series of papers are supported by 
the obtained differences between the electron temperature 
gradients of PN and HII regions, in view of the fact that 
these gradients essentially reflect the radial abundance gra- 
dients in the galactic disk. In order to make a direct com- 
parison with the estimated flattening rate of the abundance 
gradients derived by Maciel et al. (2005J) , we can convert 
the T e gradient into the equivalent O/H gradient. For HII 
regions we can use the calibration by Shaver et al. 0-983), 
according to which the oxygen gradient is related to the 
electron temperature gradient by 



dlog(Q/H) 
dR 



-1.49 x 10" 



dR 



(1) 



where the temperature gradient is in K/kpc and the oxy- 
gen gradient is in dex/kpe. This relation is also supported 
by more recent work on HII regions, such as the analysis 
by Deharveng et al. (2000). The oxygen gradient for HII 
regions is then d\og(0/H)/dR ~ —0.043 dex/kpc, which is 
similar to the value obtained by Deharveng et al. (2000) 
based on an entirely different sample. Also, an O/H gradi- 
ent of —0.044 dex/kpc was recently obtained by Esteban 
et al. (|2005p from oxygen recombination lines, a method 
almost independent of the assumed electron temperatures, 
and totally independent of the [OIII] forbidden lines. A sim- 
ilar estimate has been presented by Quireza et al. ([2006), 
also based on the Shaver et al. (1983) calibration. 

For planetary nebulae, we can have an idea of the cor- 
responding O/H gradient by inspecting Table 1 of Maciel 
et al. (|2005p for Group II objects (ages of 4-5 Gyr), from 
which we get dlog(0/H)/dR ~ -0.089 ± 0.003 dex/kpc. 
Alternatively, we can compute the O/H gradient directly 
for the PN sample adopted here, in which case we get a simi- 
lar value, dlog(0/H)/dR ~ -0.090 dex/kpc. Therefore, the 
flattening rate of the oxygen gradient can be estimated by 
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so that \ — 0.0094 dexkpc -1 Gyr -1 where we have used 
At ~ 5 Gyr. In view of our discussion on the electron tem- 
peratures of planetary nebulae, we would expect the rate 
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to be somewhat smaller than this value, which is in excel- 
lent agreement with the range estimated by Maciel et al. 
(|2005|) . namely x - 0-005 - 0.010 dexkpc -1 Gyr -1 . 

As mentioned in the introduction, abundance gradients 
and their time variation are valuable constraints for chem- 
ical evolution models. As an illustration, we have com- 
pared our derived flattening rate with the predictions of 
some recently published models for the Milky Way. As dis- 
cussed by Maciel et al. (|2006p there may be large discrep- 
ancies between different chemical evolution models, even 
whithin the so-called 'inside-out' class of models. In partic- 
ular, Hou et al. ([2000) adopted an exponentially decreasing 
infall rate for the galactic disk, in which a rapid increase 
in the metal abundance at early times in the inner disk 
leads to a steep gradient. As times goes on, the star for- 
mation migrates to the outer disk and metal abundances 
are enhanced in that region, with the consequence that the 
gradients flatten out. A rough estimate for the O/H gradi- 
ent variation in these models leads to a steepening rate of 
X 0.0040 - 0.0060 dexkpc -1 Gyr -1 , which is consistent 
with our present results. A similar behaviour has also been 
obtained by Alibes et al. (|2001|) . On the other hand, models 
such as those based on two infall episodes by Chiappini et 
al. (|200ip . lead to some steepening of the gradients, even 
though the inside-out approach is adopted. Possibly, the 
main reason for the different predictions of the quoted mod- 
els appears to reside on the different adopted timescales for 
star formation and infall rate, so that we expect our present 
results may be helpful in order to constrain these quantities. 
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